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ABSTRACT 

We present Chandra and XMM-Newton X-ray data of NGC 5253, a local starburst- 
ing dwarf elliptical galaxy, in the early stages of a starburst episode. Contributions to 
the X-ray emission come from discrete point sources and extended diffuse emission, 
in the form of what appear to be multiple superbubbles, and smaller bubbles prob- 
ably associated with individual star clusters. Chandra detects 17 sources within the 
optical extent of NGC 5253 down to a completeness level corresponding to a luminos- 
ity of 1.5 X 10"^^ erg s~^ . The slope of the point source X-ray luminosity function is 
— 0.54±q'^|, similar to that of other nearby dwarf starburst galaxies. Several differ- 
ent types of source are detected within the galaxy, including X-ray binaries and the 
emission associated with star-clusters. 

Comparison of the diffuse X-ray emission with the observed Ha emission shows 
similarities in their extent. The best spectral fit to the diffuse emission is obtained with 
an absorbed, two temperature model giving temperatures for the two gas components 
of ~ 0.24 keV and ^ 0.75 keV . The derived parameters of the diffuse X-ray emitting 
gas are as follows: a total mass of ~ 1.4 x 10^/^^^ Mq , where / is the volume filling 
factor of the X-ray emitting gas, and a total thermal energy content for the hot X- 
ray emitting gas of 3.4 x 10^''/-'^/^ erg. The pressure in the diffuse gas is P/k ~ 
10^/~^/^ K cm~'^. We find that these values are broadly commensurate with the mass 
and energy injection from the starburst population. Analysis of the kinematics of the 
starburst region suggest that the stellar ejecta contained within it can escape the 
gravitational potential well of the galaxy, and pollute the surrounding IGM. 

Key words: ISM: jets and outflows - galaxies: individual: NGC 5253 - galaxies: 
starburst - X-rays: galaxies. 



1 INTRODUCTION 

NGC 5253 is a starbursting dwarf elliptical galaxy which lies 
at a distance of ~ 3.15 Mpc (Freedman et al. 2001) and 
with an inclination of 67° (the mean value quoted by the 
Lyon/Meudon Extra-galactic Database [LEDA]). The signa- 
ture of Wolf-Rayet stars have been detected in the nucleus 
of the galaxy (Schaerer et al. 1997, Walsh & Roy 1987) im- 
plying that this is a young starburst and as such allows an 
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opportunity for the study of the starburst phenomenon in its 
earlier stages of development and its effect on galaxy evolu- 
tion in the local Universe. In fact, Rieke, Lebofsky & Walker 
(1988) went so far as to classify it as one of the youngest star- 
bursts known. Dwarf galaxies, as the basic building blocks 
in the hierarchical merging cosmology scenario, are likely 
to have harboured the earliest sites of star-formation in the 
Universe and so their study in the local Universe can give in- 
sight into the evolution of such objects at high redshift. Ob- 
servations of local edge-on starburst galaxies (Strickland et 
al. 2000; Weaver 2001) are presenting a picture of kpc-scale, 
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soft X-ray emitting, bipolar outflows in the form of galactic 
winds transporting mass, newly synthesised heavy elements 
and energy into the intergalactic medium (IGM). A simi- 
lar outflow is seen in the dwarf starburst galaxy NGC 1569 
(Martin, Kobulnicky & Hcckman 2002), while the dwarf 
NGC 4449 shows what may be the beginnings of a galac- 
tic wind emerging from an extended superbubble (Summers 
et al. 2003). The same situation seems to be the case in 
the NGC 3077 dwarf (Ott, Martin & Walter 2003). These 
winds result from the pressure driven outflows along these 
galaxy's minor axis produced from the efficient thermaliza- 
tion of the mechanical energy from the supernovae (SN) ex- 
plosions and stellar winds of the massive stars in their OB 
associations and super star-clusters (SSC). NGC 5253 being 
both inclined to our line-of-sight and young presents a less 
clearly observable picture. Its youth means that it is un- 
likely to have developed a superwind and the diffuse X-ray 
emission observed is much more likely to be associated with 
multiple superbubbles around its OB associations and SSC, 
as suggested from earlier i^OiSAT observations of this galax:y 
(Strickland & Stevens 1999). 

Observations of NGC 5253 at radio wavelengths have 
shown it to have a very flat centimetre-wavelength contin- 
uum which is indicative of thermal emission from HII re- 
gions, with only low levels of synchrotron emission from su- 
pernova remnants (Beck et al. 1996; Turner, Ho & Beck 
1998). Ferrarese et al. (2000) quote a low metallicity of 
0.2^0, based on the work of Webster & Smith (1983). This 
figure suggests that a large amount of metal enrichment 
has not occurred in the central HII regions, and may be 
further evidence of the youth of the starburst region. The 
total and HI masses for the galaxy are calculated to be 
6.4 X 10* Mq and 8.3 x lO'' M© from 21 cm observations 
(Reif et al. 1982), corrected for the difference in assumed 
distance to NGC 5253. Later VLA 21 cm observations show 
NGC 5253 to be peculiar, in so much as its neutral hydrogen 
appears to rotate about the optical major axis of the galaxy 
(Kobulnicky & Skillman 1995). However, CO observations 
(Turner, Beck & Hurt 1997; Meier, Turner & Beck 2002) 
detect sources which are coincident with the dust lane seen 
to the SE of the nucleus of the galaxy in optical images, and 
these CO clouds appear to be infalling into NGC 5253. Con- 
sequently, the dynamical situation in NGC 5253, rotation or 
infall, is not completely clear. 

Infrared emission is detected from what seems to be a 
highly obscured, massive (10-^ - 10*^ stars), small (~ a few 
pc diameter) and very young globular cluster in the cen- 
tral starburst region of the galaxy (Gorjian, Turner & Beck 
2001). This source contributes ~ 50% of the total observed 
infrared luminosity of the galaxy. A summary of some of 
the properties of NGC 5253 is given in Table 1, along with 
estimates of the current star-formation rate of the galaxy 
estimated from far-infrared and Ha fluxes. 

Caldwell & Phillips (1989), using broadband multi- 
colour (UBV) CCD images, narrow-band Ha and long-slit 
Call triplet spectra, determined that NGC 5253 had un- 
dergone an increased rate of star-formation throughout the 
galaxy about 10* — 10® yr ago, as suggested by the presence 
of over 100 star clusters of a similar age lying outside the 
nuclear region. They also established that the intense star- 
formation now occurring in the nuclear region of the galaxy 
began around 10 Myr ago, giving an upper limit for the age 



of superbubbles associated with it. More recently, Tremonti 
et al. (2001) analysed several clusters in NGC 5253, deter- 
mining ages in the range of 1 to 8 Myr. Calzetti et al. (1997) 
also noted two other clusters with older ages, determined to 
be between 10 and 50 Myr. 

The two very bright SN explosions seen in this galaxy 
in 1895 and 1972 are amongst the apparently brightest 
recorded extrargalactic SNe (Ardeberg & de Groot 1973) 
and may be associated with the large scale star-formation, 
although neither was particularly centrally located. The 
Fabry-Perot Ha images of Marlowe et al. (1995) show that 
the ionized gas in the galaxy has a complex distribution, 
consisting of a bright central region extending ^ 30" from 
the centre in a N-S direction, embedded within a system of 
loops and filaments, with an overall extent of ~ 2'. To the 
NW and WSW of the centre are 2 large superbubbles with 
diameters ~ 1', expanding at ~ 35 km s~^ , while to the east 
side of the galaxy, a pair of radially oriented filaments are 
detected. This latter observation is suggestive of the rupture 
of a superbubble blown by the stellar activity in the central 
regions. 

NGC 5253 has previously been observed at X-ray wave- 
lengths with the Einstein IPC (Fabbiano, Kim & Trincliicri 
1992) and both the ROSAT PSPC (Martin & Kcnnicutt 
1995; Stevens & Strickland 1998) and HRI (Strickland & 
Stevens 1999). The Einstein IPC and ROSAT PSPC results 
showed soft thermal X-ray emission from what appeared to 
be an extended source, suggestive of a luminous superbub- 
ble. The better spatial resolution of the ROSAT HRI instru- 
ment detected a complex of at least five X-ray sources, that 
could be young superbubbles blown by individual young star 
clusters in the starburst region. The total X-ray luminosity 
determined from these observations was similar and lies in 
the range of (2.4 — 4.1) x 10^* erg s"'^ , depending on the 
value adopted for the absorbing column density. 

In Section 2 we describe the Chandra and XMM-Newton 
observations. The X-ray emission from the whole galaxy, 
the point sources and the diffuse emission axe discussed in 
Section 3. Section 4 contains a more general discussion of 
the morphology of the X-ray emission and its relationship to 
emission from other wavebands, along with a determination 
of the dynamics of the expanding superbubbles of the central 
region and their effect on NGC 5253. Our main conclusions 
are summarised in Section 5. 



2 OBSERVATIONS AND ANALYSIS 

A 57 ks Chandra observation of NGC 5253 was obtained on 
Jan 13*'' - 14*'' 2001 and a 46 ks XMM-Newton observa- 
tion was obtained on 8*'' — 9*'' August of the same year. 
For the Chandra data, analysis was carried out on the data 
contained within the S3 chip of the ACIS— S instrument us- 
ing CIAO (version 2.2.1), HEASOFT (version 5.1), XSPEC 
(version 11.1.0) and ASTERIX (version 2.3-bl). The data 
were first reprocessed using CALDB version 2.9 and then 
filtered to remove periods of flaring and lower than aver- 
age count rates using the lc_clean.sl script with a 5a clip- 
ping, leaving a total of 47.8 ks of useful data. The data were 
then further filtered to retain only data in the energy band 
from 0.3 — 8.0 keV and were background subtracted using 
the appropriately scaled Chandra X-ray Observatory Center 




Figure 1. Low resolution (smootlied using a Gaussian with FWHM of 4 pixels ~ 2 ), background subtracted Chandra image, in the 
0.3 — 8.0 keV energy band, showing the full ACIS-S3 chip field of view, marked with the 31 point sources (Table 2). The image is an 
inverted log greyscale with the flux density ranging from 7.63 X 10~^* erg cm~^ arcmin"^ to 9.72 X 10~^^ erg s~^ cm"'^ arcmin"'^. 
The D25 ellipse is also shown. North is to the top and East is to the left. 



(CXC) background event data set. After production of an 
exposure map, the point sources present in the data were 
detected using the CIAO wavdetect tool, run with an expo- 
sure map and a source significance threshold of 9.5 x 10~^, 
a value of ~ (number of pixels) which should have limited 
the number of false detections to ^ 1. A total of 31 point 
sources were detected and these are shown in Fig.0 overlaid 
on a smoothed image of the ACIS-S3 chip (smoothed with a 
Gaussian having a FWHM of 4 pixels ~ 2"), while their po- 
sitions and background subtracted count rates are listed in 
Table 2. Of these 31 sources, 17 lie within the optical extent 
of the galaxy, as shown by the D25 ellipse (de Vaucouleurs 
et al. 1991) on Fig. [T] These 17 sources are the ones most 
likely to be associated with NGC 5253 and so are the only 
ones considered further in this analysis. We note that results 
from this data have already been presented as part of a sur- 
vey of the point source populations of nearby galaxies by 
Colbert et al. (2004), and a comparison with these results 
will be presented later. 

From the XMM-Newton observation, only the European 
Photon Imaging Camera (EPIC) data has been analysed. 
Removal of flares using a recursive 3a clipping technique 



and standard pattern and flag flltering for the three cam- 
eras left 39.0 ks of useful data for MOSl, 39.6 ks for M0S2 
and 26.7 ks for PN. These data were then also filtered to 
retain only data in the 0.3 — 8.0 keV energy band in line 
with the Chandra data. A background subtracted (using 
the David Lumb 'blank sky' and Philippe Marty 'closed' 
events files) mosaic of the exposure corrected images from 
the three cameras is shown in Fig. |5| The (0,0) point on 
this image lies at the centre of the D25 ellipse at a position 
of a = 13''39"'56.2", 5 = -31°38'29.9". The diffuse X-ray 
emission is seen to be extended to the WSW of the nucleus 
which is not apparent in the Chandra image. The latter does 
however show evidence of absorption of the X-ray emission 
due to the presence of the dust lane that extends E from the 
nuclear region, which is not apparent on the XMM-Newton 
data. It is also apparent from comparing Figs. Qand|21that 
some of the sources within the Chandra ACIS-S3 field-of- 
view are exhibiting variability. In the eight months between 
the two observations, sources 25 and 30 have faded whilst 
six possible additional sources are apparent in the XMM- 
Newton data, to the east of source 1, south-east of source 7, 
north-east of source 12, north of source 24, north of source 
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Figure 2. Low resolution (smoothed using a Gaussian with FWHM of 2 pixels '^8 ), XMM-Newton, background subtracted and exposure 
corrected, mosaiced image of the data from the EPIC MOSl, MOS2 and PN cameras, in the 0.3 — 8.0 keV energy band. The image is an 
inverted log greyscale with the flux density ranging from 5.2 X 10~^*erg s~^ cm~^ arcmin"^ to 2.1 X lO^^'^erg s~^ cm~^ arcmin"^, 
with North to the top and East to the left. The D25 ellipse is also shown. The X-ray emission is extended towards the WSW. 



26 and north of source 27. The last three in this list may 
be associated with NGC 5253 as they lie just inside the D25 
ellipse. Fig. 13 shows the position and extent of the Chandra 
sources overlaid on the unsmoothed XMM-Newton mosaic 
and shows that a large fraction of what appears to be ex- 
tended diffuse emission to the SW along the galaxy's major 
axis is more likely to be associated with point sources. How- 
ever the diffuse emission does show extension to both the 
NW and SE along the minor axis, which may be associated 
with the shells and radial filaments seen in Ha images (Mar- 
lowe et al. 1995). 

After detection, the point sources were subtracted from 
the Chandra data so that the diffuse X-ray emission con- 
tained within that observation could be investigated. Also, 
a 3 colour adaptively smoothed image of the central region 
of the galaxy has been produced using csmooth in fft (fast 
Fourier transform) mode and this is shown in Fig. |1] The 
lower and upper sigmas for deriving the smoothing kernel 
were set at 2 and 5 respectively. This image shows that the 
spectral characteristics of the sources vary and that the dif- 
fuse emission appears to consist mainly of soft emission lying 
behind regions if varying absorption. The red, green and blue 



images used correspond to energy bands of 0.3 — 1.1 keV , 
1.1 — 3.8 keV and 3.8 — 8.0 keV respectively. These bands 
were chosen so that the total number of counts in each band 
was approximately equal and thereby maximising the S/N 
ratio for the bands. The diffuse emission is seen to be con- 
fined predominantly to the nuclear region of the galaxy and 
the position of the indentation seen to the East, suggesting 
heavy absorption, corresponds to the position of the dust 
lane identified in the galaxy (Hunter 1982). In addition the 
effect of the dust is seen, in the reduction of intensity and 
change of colour as a result of the absorption of the softer 
emission, in the darker band running from E to W across 
the nuclear region. 

3 RESULTS 

3.1 Integrated Emission 

The background subtracted spectrum of the total X-ray 
emission from within the D25 ellipse of NGC 5253 was fitted 
separately for the four instruments (ACIS-S, MOSl, M0S2 
and PN), using the modified Levenberg-Marquardt method 
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Table 1. Details of NGC 5253. The distance is from Freedman 
et al. (2001), other data is from the LEDA and NED databases. 
The infrared luminosity Ljr is calculated using the 12, 25, 60 
and lOO^tm IRAS fluxes (Sanders & Mirabel 1996) and the Ha 
luminosity is from Marlowe et al. (1997, corrected for distance). 
The corresponding star-formation rates (SFR) for both the IR 
and Ha luminosities are calculated using the conversion formulae 
in Kennicutt (1998). 



Parameter 


Value 


Classification 


Im pec, HII 


Distance 


3.15Mpc 


RA (J2000) 


12 39 56.0 


Dec (J2000) 


-31 38 36 


D25 ellipse 


4.8' X 1.9' 


Position angle 


45° 


niB 


10.75 


B-V 


0.43 




4.5 X 10*2 erg s"! 




2.6 X 10*° erg s"! 


SFR rate (IR) 


0.2 M0 yr-l 


SFR rate (Ha) 


O.2M0 yr-l 



Table 2. Positions and count rates of the 31 sources detected 
in the NGC 5253 Chandra S3 chip data. Column 1 is the source 
number ordered in increasing R.A.. Columns 2 and 3 give the 
R.A. and Dec. of each source and Column 4 lists their background 
subtracted count rates. 



Source RA (h m s) Dec (° ' ") Count Rate 

(xlO-3 cts s^i) 



1 


13 


39 


37.73 


-31 


40 


15.8 


0.77 ±0.15 


2 


13 


39 


39.16 


-31 


36 


06.9 


0.32 ±0.09 


3 


13 


39 


39.17 


-31 


36 


09.2 


0.42 ±0.10 


4 


13 


39 


43.15 


-31 


39 


20.2 


1.07 ±0.16 


5 


13 


39 


45.12 


-31 


36 


58.2 


0.37 ±0.10 


6 


13 


39 


45.53 


-31 


42 


06.1 


0.39 ±0.10 


7 


13 


39 


46.69 


-31 


42 


28.0 


0.43 ±0.11 


8 


13 


39 


47.27 


-31 


41 


40.7 


1.95 ±0.22 


9 


13 


39 


50.66 


-31 


39 


20.7 


1.50 ±0.18 


10 


13 


39 


51.17 


-31 


38 


40.8 


0.32 ±0.09 


11 


13 


39 


51.62 


-31 


40 


16.2 


0.48 ±0.11 


12 


13 


39 


51.96 


-31 


41 


43.9 


4.03 ±0.30 


13 


13 


39 


52.07 


-31 


38 


59.7 


0.19 ±0.07 


14 


13 


39 


53.36 


-31 


43 


43.9 


1.42 ± 0.20 


15 


13 


39 


53.80 


-31 


38 


32.7 


0.66 ±0.13 


16 


13 


39 


55.39 


-31 


38 


32.8 


0.80 ±0.17 


17 


13 


39 


55.69 


-31 


38 


31.5 


1.02 ±0.20 


18 


13 


39 


55.85 


-31 


38 


22.8 


1.03 ±0.20 


19 


13 


39 


56.04 


-31 


38 


24.8 


1.43 ±0.23 


20 


13 


39 


56.15 


-31 


38 


35.9 


0.58 ±0.16 


21 


13 


39 


56.35 


-31 


38 


20.5 


8.22 ±0.43 


22 


13 


39 


56.45 


-31 


38 


25.8 


1.53 ±0.21 


23 


13 


39 


56.67 


-31 


42 


45.3 


1.57 ±0.21 


24 


13 


39 


57.43 


-31 


37 


32.5 


1.50 ±0.06 


25 


13 


39 


59.27 


-31 


37 


50.2 


0.33 ± 0.09 


26 


13 


40 


02.19 


-31 


37 


19.1 


0.96 ±0.15 


27 


13 


40 


03.40 


-31 


37 


14.5 


0.78 ±0.13 


28 


13 


40 


03.84 


-31 


44 


03.7 


1.25 ±0.21 


29 


13 


40 


04.29 


-31 


37 


30.3 


0.41 ±0.10 


30 


13 


40 


04.97 


-31 


42 


03.0 


0.24 ±0.08 


31 


13 


40 


10.81 


-31 


40 


33.2 


0.29 ±0.09 




Figure 3. Unsmoothed background subtracted and exposure cor- 
rected, mosaiced image of the data from the XMM-Newton EPIC 
MOSl, MOS2 and PN cameras in the 0.3-8.0 keV energy range. 
The overlaid small ellipses show the positions and sizes of the 
sources detected in the Chandra data. The D25 ellipse is also 
shown. The extension in the X-ray emission towards the WSW 
is largely associated with point sources, but the emission is also 
slightly extended to the NW and SE along the galaxy's minor 
axis, as well as along the major axis. North is to the top and East 
is to the left. 

and standard statistic from XSPEC, with an absorbed 
2 thermal component plus power-law fit, the same model 
that was fitted in earlier work to the similar dwarf star- 
burst galaxy NGC 4449 (Summers et al. 2003) . It was found 
though, for NGC 5253 that each component had to be fitted 
with a separate absorbing column, in addition to the absorp- 
tion due to the interstellar medium within the Milky Way 
(Galactic column density, Nh(gal) = 3.87 x 10^" cm~^), 
to obtain a satisfactory fit. Another difference in the fitting 
for these two galaxies was in the energy range over which 
the spectral fitting was carried out. The NGC 5253 Chan- 
dra data was found to contain both a hard excess above 
~ 6.0 keV and a soft excess below ~ 0.37 keV , neither of 
which could be adequately dealt with by either the flare re- 
moval or background subtraction. Even the use of a local 
background instead of the background event data set did 
not remove them and so the spectral fitting for all the in- 
struments was confined to the 0.37 — 6.0 keV energy band 
for ease of comparison of the results from the two differ- 
ent telescopes. In addition, a single simultaneous fit to the 
data from the three XMM-Newton cameras was performed. 
The two absorbed thermal components were modelled using 
the wabs (photo-electric absorption using Wisconsin cross- 
sections) and mekal thermal plasma codes within XSPEC. 
Initially, Galactic absorption was assumed for all the column 
densities and a value of 0.19 Zq assumed for the metallic- 
ity of the X-ray emitting gas (Ferrarese et al. 2000). The 
results are summarised in Table 3 and the fitted spectra 
for the Chandra data and simultaneous fit to the XMM- 
Newton data are shown in Fig.|3 Two thermal plasmas were 
used to represent the diffuse emission firstly because multi- 
phase models are seen to be needed to best fit the emission 
from other starburst galaxies (e.g. NGC 253, Strickland et 
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al. 2002; NGC 1569, Martin et al. 2002; NGC4449, Summers 
et al. 2003), and secondly because single temperature fits to 
just the diffuse emission do not have a wide enough energy 
distribution to fit the spread seen in the data and so are 
statistically less robust. 

We note that the MOSl data gives substantially differ- 
ent results from the other two XMM-Newton cameras, par- 
ticularly for the column densities and metallicities. There 
is also a difference between the results obtained from the 
two satellites, with the Chandra data giving higher col- 
umn densities for the soft thermal and for the power-law 
components. The average unabsorbed flux, in the 0.37 — 
6.0 keV energy band, from NGC 5253 for the Chandra data 
is (1.56±o;i2) X 10~^^erg s~^ cm~^, and for the combined 
XMM-Newton data is (3.22±!]:|o) x 10'" erg s~^ cm~^. 
These values correspond to a total intrinsic X-ray luminos- 
ity of (1.86±o:") X 10^^ erg s'^ from the Chandra data and 
(3.83±oi7) X 10^^ erg s"^ from the XMM-Newton data for 
our assumed distance of 3.15 Mpc. In both cases, these val- 
ues are higher than the results obtained from earlier ROSAT 
and Einstein data but the differences are explained by the 
higher fitted column densities of the more recent data. The 
fraction of the emission due to the individual components is 
15% - cooler thermal, 25% - hotter thermal and 60% - power- 
law from the Chandra data and 40% - cooler thermal, 50% - 
hotter thermal and 10% - power-law from the XMM-Newton 
data. The very high fitted column density for the power-law 
component in the Chandra data accounts for a lot of this 
discrepancy. The differences in the calibration between the 
MOSl camera and the other two XMM-Newton cameras will 
also lead to possible discrepancies when comparing results 
from the combined XMM-Newton fits to Chandra fits. We 
do however note the Chandra calibration problems at low 
energy, which may affect matters at low energies. An impor- 
tant point is that that the discrepancies between fits to the 
XMM-Newton and Chandra data are much smaller for the 
two thermal components, and it is these components that 
we concentrate on in this paper. 

3.2 Point Sources 

A background subtracted spectrum was extracted from the 
Chandra data for each of the 17 sources within the D25 el- 
lipse. (The background spectra were taken from the afore- 
mentioned CXC background file and were appropriately 
scaled. The use of the background files meant that the back- 
ground spectrum was taken from the identical position on 
the CCD and avoided the problem of contamination from 
adjacent sources, which as shown in Fig. is particularly 
a problem in the central region of NGC 5253.) The spec- 
tra were grouped so that each had a minimum of 5 data 
bins and when possible a minimum of 10 counts per bin. 
(In fact, due to the low statistics, this was only possible 
in the case of sources 9, 19, 21 and 22.) No attempt was 
made initially to fit the data where a source had less than 
50 counts after background subtraction. Source 21 has the 
highest count-rate by far of any of the sources (see Table 2) 
and is fitted well by an absorbed two component (mekal 
plus power-law) fit. Galactic absorption is assumed to be 
Nh(gal) = 3.87 X 10^° cm"^ and the metallicity of the 
X-ray emitting gas is fixed at O.13Z0 [the average value fit- 
ted to the total X-ray emission from NGC 5253, which is 



slightly lower than, but agrees within errors with, the value 
of 0.19 ±Si ;o| Zq obtained for the HII regions in the cen- 
tre of NGC 5253 by Webster & Smith (1983)]. The fitted 
spectrum for this source is shown in Fig. |S| and the fitted 
parameters are summarised in Table 4 along with those for 
the other 16 sources within the D25 ellipse. Where possible, 
the other sources were fitted with this same model unless 
the normalisation went to zero for one of the components, 
in which case that component was omitted and these fits are 
indicated by a dash in the parameter values in Table 4. As 
suggested by Strickland & Stevens (1999), from analysis of 
the ROSAT HRI data for this galaxy, some of the sources 
detected in the central region of NGC 5253 may be due to 
superbubbles surrounding super star-clusters rather than in- 
dividual point sources, such as X-ray binaries. Two of these 
five sources found by Strickland & Stevens (1999) have been 
further resolved into two sources using the superior resolu- 
tion of Chandra. We can relate the HRI sources recorded in 
Strickland & Stevens (1999) to the Chandra sources found 
here as follows: fJOSAT source A corresponds to source 21, 
B to source 22, C to source 20, D to sources 18/19 and 
E to sources 16/17. The assertion that these sources could 
be extended superbubbles is supported by the source detec- 
tions here where in particular, sources 16, 17, 18, 19 and 20 
are all detected as being elliptical whereas the PSFs at the 
corresponding positions on the CCD are all circular. The in- 
clusion of a thermal component in the spectral fits for these 
sources is then to be expected and in the cases of sources 
16, 17, 18 and 20, the thermal components completely domi- 
nate the spectra. Source 21 appears to be substantially more 
luminous in this Chandra observation than in the ROSAT 
HRI observation, though a detailed comparison is difficult 
because of the low source counts in the HRI observation. 

From the best fits to all the sources within the D25 el- 
lipse, the absorption corrected fluxes and luminosities in the 
0.3 — 8.0 keV energy band were determined. These luminosi- 
ties are also included in Table 4 and the log(Ai') — \og{Lx) 
plot for the 17 sources possibly associated with NGC 5253 is 
shown in Fig.|7| This has been fitted with a single power-law, 
with a slope of — 0.54±o;ig for the higher luminosity sources. 
Sources with luminosities below an absorption corrected lu- 
minosity of 1.46 X 10'^^ erg s~^ were not included in the fit, as 
these figures represent the completeness limit for the data. 
The sources where the thermal component dominates are 
included in this data, even though these may represent the 
emission from young superbubbles rather than point sources. 
The slope determined for all the sources is not quite as steep 
(though comparable to within the errors) as the —0.70 found 
for NGC 4214 (HartweU et al. 2003) but slightly steeper than 
the -0.51 found for NGC 4449 (Summers et al. 2003). As a 
further comparison, Martin et al. (2002) give an upper limit 
of —0.55 for the slope of the luminosity function of sources 
detected in the 1.1 — 6.0 keV energy band in NGC 1569, but 
also state that the value is not well constrained. The low 
source counts in the NGC 5253 data make anything more 
rigorous than a statement of the slope values impossible. The 
NGC 5253 value is steeper than the value quoted (~ —0.45) 
by Zezas et al. (2001) for the larger starburst galaxies, M82 
and the Antennae. The steeper slopes seen in the dwarfs may 
be the result of them having fewer high luminosity sources, 
and in the case of NGC 5253 this could be directly related 
to it being a younger starburst. (For more analysis of the 
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Figure 4. Adaptively smoothed 3 colour image of a 3 arcmin square region of the Chandra ACIS-S3 chip, centred on the nucleus of 
NGC 5253 at a = 13''39"'56.2^ <5 = -31°38'29.9" (red: 0.3-1.10 koV , green: 1.10-3.80 keV and blue: 3.80-8.0 keV). The hardness 
variation in the source spectra is evident, while the diffuse emission is seen to be predominantly soft. 



Table 3. Best fits to the background subtracted spectra of the total X-ray emission from within the D25 ellipse of NGC 5253 for a 
wabsQjXi^ {wabs{mekal) + wabs{mekal) + wabs{po)) model, with the Galactic column wabsQ^i^ = 3.87 X 10^" cm~^. The five columns of 
results proceeding from left to right are for the Chandra ACIS-S3 chip, the individual XMM-Newton MOSl, MOS2 and PN cameras and 
the combined results from the data for the three XMM-Newton cameras fitted simultaneously. Row 1 lists the individual instruments, 
rows 2, 4, and 7 are the fitted column densities for the three components, Af^f^jj, Nu(2) ^i^d ^H{3) respectively. Rows 3 and 5 are 
the temperatures, kTi and kT2 of the two thermal components and row 6 is the fitted metallicity of the X-ray emitting gas, fitted 
simultaneously for the two thermal components. Row 8 gives the photon index, F, of the power-law component, while row 9 gives the fit 
statistics. 



Chandra XMM-Newton 
Instrument ACIS-S MOSl MOS2 PN 3 XMM 



A^(i) (xl02i cm-2) 

fcTi ( keV ) 

Nh(2) (xl02i cm-2) 

kT2 ( keV ) 

Z (Zq) 

Nh(3) (xl02i cm-2) 

r 

xVd.o.f. 



1.01±0:J2 l.00±0;i 0.89±0:3i 0.79±0:12 o.48±0:09 

0.25 ±0.01 0.26 ±0.02 0.22±g gf 0.22 ± 0.01 0.25 ± 0.01 

3.84±0:34 2.06±0:|I 1.34±o;jo 0.82±0:2| 3.78±0:i| 

0.76 ±0.04 0.64 ±0.08 0.59±j;-g|j 0.73±g [|^ 0.67±{J g| 

0.13 ±0.01 0.16 ±0.02 0.11 ±0.01 0.13 ±0.01 0.11 ±0.01 

46.4±|:M 0.51±J:|t 5.00±»giO 8.77±I:?1 7.01±i«32« 

1.69±0:06 3.60±;;IB 2.12±i:6i 1.92±0:g 1.54±0:i2 

206/150 66/53 71/66 159/145 323/284 



luminosity functions of these and other starburst galaxies 
see Hartwell et al. 2003). The efltect of omitting the sources 
dominated by thermal emission which are probable young 
superbubbles (sources 16, 17, 18, 19 and 20) is to decrease 
the slope of the fit to —0.41, which is more in line with 
the other results quoted above. This highlights the fact that 
avoiding source confusion is crucial to this type of analysis 
and the inclusion of super star-clusters within superbubbles. 



as point sources biases the results in such a way as to steepen 
the slope of the point source luminosity function. 

We note that Colbert et al. (2004) in their analysis 
of the NGC 5253 data, quote a steeper value for the point 
source luminosity function of —0.92. A major difference be- 
tween this and that analysis is that here we fit the spectrum 
of each source individually (Table 4), whereas Colbert et al. 
(2004) determine source fluxes using a power-law model with 
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Figure 5. Fitted spectrum of the total X-ray emission from NGC 5253 Left panel: Chandra data, right panel: simultaneous fit to the 
data from the three XMM-Newton cameras. The fit is for an absorbed two temperature plus power-law fit, with a column density of 
^H(GAL) — 3.87 X 10^'^ cm~^ assumed for Galactic absorption. See Table 3 and text for details of the individual fits. The individual 
contributions from the separate components are shown on the Chandra data. 



a fixed slope and a column set at the Galactic absorption 
level. This may well explain the different slopes determined 
for NGC 5253. 

The data were point source searched in the three dif- 
ferent energy bands used to produce the three colour im- 
age, (soft 0.3 — 1.1 keV , medium 1.1 — 3.8 IceV and hard 
3.8 — 8.0 keV ) and in the three energy bands used for 
the analysis of NGC 4449 (soft 0.3 - 0.8 keV , medium 
0.8 - 2.0 keV and hard 2.0 - 8.0 keV ) so that a direct com- 
parison could be made of the sources in the two galaxies. 
The hardness ratios for the individual sources within the 
7325 ellipse were calculated for the sources detected in more 
than one of the energy bands. Due to the low count rates 
for many of the sources most were detected in only 1 or 2 
bands. The total counts in each of these energy bands and 
the calculated hardness ratios are shown in Tables 5 and 6, 
for the two different sets of energy bands. Note, that when 
point source searching was carried out in the separate en- 
ergy bands, sources 18 and 19 when detected as a single 
source rather than two sources as in the broadband search- 
ing. The two hardness ratios shown are (m — s)/ (m + s) and 
{h — m)/{h + m), where s, m and h are the counts in the 
soft, medium and hard bands respectively. It is noticeable, 
particularly in Table 5, that the sources in the central re- 
gion (sources 16-20) are all very soft, with very few counts 
above 1.0 keVand virtually none above 2.0 keV , suggesting 
that these are associated with the emission from superbub- 
bles surrounding super star-clusters rather than individual 
harder objects like X-ray binaries. 

In general, the sources detected in NGC 5253 are softer 
than those in NGC 4449 and lie behind an on average higher 
absorbing column density. The implication of this being that 
the sources in NGC 5253 must be softer as the higher col- 
umn density would preferentially absorb the softer X-ray 
emission. The low counts for the sources in NGC 5253 will 
effect these results too, but still, only 5 out of 17 sources 
have detections in the 2.0 — 8.0 keV band compared with 
15 out of 24 for NGC 4449. To better compare some of these 
sources, the sources where hardness ratios can be calculated 
from the NGC 5253 data, in the 0.3 - 0.8 keV ,0.8-2.0 keV 




1 2 5 

channel energy (keV) 



Figure 6. Spectral fit obtained for source 21, the most luminous 
source in the Chandra data. The fit is an absorbed two component 
(mekal+power-law) fit with the Galactic column density assumed 
to be 3.87 X 10^° cm~2 and the metallicity fixed at 0.13^0 (the 
average fitted value from the fits to the total emission). The best 
fit gives Nh = (8.56±[;-|^) X lO^l cm-^, kT = 0.24 ± 0.01 keV 
and the photon index, f = 2.09±° }^. 



and 2.0 — 8.0 keV energy bands, have been overlaid on the 
hard vs. soft hardness ratio plot for NGC 4449. The results 
are shown in Fig. |S| where the NGC 4449 data is shown in 
black and that for NGC 5253 is in red. The sources repre- 
sented by squares with arrows, to the right and bottom of 
the plot, were only detected in 2 of the 3 energy bands and 
so their positions represent the most extreme positions pos- 
sible for them on this plot. 

The average fitted column density for the NGC 4449 
sources is ~ 3.0 x 10^^ cm"^ while that for the NGC 5253 
sources is ~ 7.5 x 10^^ cm~^. The effect of these differing 
values on a plot such as Fig. |H| is to displace the sources 
towards the right and slightly towards the top of the plot, 
as the absorbing column behind which they lie increases. 
For a truer comparison, the NGC 5253 plotted values would 
decrease by ~ 0.2 — 0.3 in the (m — s)/{m -\- s) ratio and 
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Table 4. Results of fitting a wabscAL (wabs(mekal+po)) model (the best fit model for source 21) to all 17 sources within the D25 ellipse. 
Column 1 gives the source numbers as shown on Fig.0 Column 2 contains the fitted column density for each source assuming a Galactic 



column density of Njj 



(GAL) 



3.87 X 10^" cm~^. Columns 3 and 4 contain the temperature of the mekal component and photon index 
of the power-law component respectively. Column 5 and 6 show the absorption corrected flux and luminosity for each source. The errors 
shown are the 90% confidence regions for each parameter. 



Source 



cm" 



kT 
( keV) 



Unabsorbed Flux 



(xlO- 



erg s 



Unabsorbed Luminosity 
(xlO^* erg s~l ) 



9 
10 
11 
13 
15 
16 
17 
18 
19 
20 
21 
22 
24 
25 
26 
27 
29 



4 4Q-I-I.24 

7.91±?6gf 

7.84±f;|i 



73 



44.59±f 
5.02±i:?4 

O.OZICq 

5.33±i-» 

6-84±[i;?i 
q qq-i-1.16 

8:ia 

'-0.50 
1 4q-|-6.42 

2.28±3-»| 
^1:33 
•=0.68 
4 5C_|_10.75 

'•-'-■^='=1.80 



8.56± 

5.35± 



0.84±A 



0.08±S 
0.25±[| 
0.23±{] 
0.25±j| 
0.25±2 
0.17±° 
0.25±[j 
0.24±2 
0.24±g 
0.22±2 
0.24±° 
0.15±° 
0.08±o 



0.18± 
0.15± 



n oq+0.07 
'^•^^='=0.07 



1.94±!J 
1 9Q-I-' 

1 45-1-0.84 

2 qq-i-1-37 



1 41+U.93 
='=0.44 

2.09±0l| 
2.50±f;g 

2.50±f|o 
0.86±i;3| 

-'-•-'-^='=0.58 

0.70±3i8 



25.7±i5.o 
0.52±o.6i 

5.50 



60 



2.10±^ 



e.ioztj 
4.oo±i:io 

4.50±HS 

9.io±i:|o 

lO.70 
10™ 

85-3±io:s 
24 q-i-6.40 

n Q7-I-O.44 
'^••^'='=0.35 
n 40-I-O.55 
'^•^•^='=0.35 

^•^'^='=1.00 

i.io±«:« 
i.6o±;:io 



n 70-1-O.66 
17-I-0" 

i-5o±g:li 

0-48±g:ii 
0.54±0:i| 

n 25-I-0 08 
in n-l-1.30 

04-1-0.05 
0.25±g:i° 

1O-I-0.05 
^•^■^=■=0.05 




Figure 7. log(A'^) — log(Lj)<; ) plot for the 17 point sources detected 
within the D25 ellipse of NGC 5253. The Lx values are absorption 
corrected and have been fitted as shown in Table 4. The line shown 
is the power-law fit to the high luminosity end of the data, with 
a slope of — 0.54±Q'^g. Incompleteness occurs at an absorption 
corrected luminosity of 1.46 X lO^'^ erg s~i . 



~ 0.1 in the [h — m)/(h + m) ratio, if the average column 
density of NGC 4449 was assumed. Of the 5 sources detected 
in the 2.0 — 8.0 keV band, two of them (sources 10 and 27) 
along with source 13, which in addition has a very high fitted 
column density, have spectra that show no emission below 
1.0 keV . As such, these three sources are possible candidates 
for background AGN. These figures are in good agreement 
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NGC 5253 red 
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Figure 8. Hard vs. soft hardness ratios for the sources detected 
in the separate energy bands. The data shown in black are for 
NGC 4449 (Summers et al. 2003) and in red for NGC 5253. The 
hard hardness ratio is calculated from the 2.0 — 8.0 keV and 
0.8 — 2.0 keV energy bands and the soft hardness ratio is cal- 
culated from the 0.8 — 2.0 keV and 0.3 — 0.8 keV energy bands, 
in the usual way, as defined in the headings of columns 5 and 6 of 
Table 6. An indication of the likely origin of the sources depend- 
ing on their location in this plot is shown by labels placed next to 
the identified sources in the NGC 4449 data. The sources plotted 
with squares to the right and bottom were only detected in 2 out 
of the 3 energy bands. The arrows indicate the direction that the 
sources would move if higher count rates permitted detections in 
all three energy bands. The error bars have also been omitted 
from these sources to avoid additional confusion on the plot. See 
Table 6 for those associated with NGC 5253. 
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Figure 9. Fitted spectrum of tlie NGC 5253 diffuse emission 
as seen with Chandra. The fit is an absorbed two tempera- 
ture plus power-law fit, with assumed Galactic column den- 
sity of A'"£f((3^£) = 3.87 X fO-^° cm~^. The fitted temperature 
components were 0.24 ± 0.01 and 0.75 ± 0.05 keV , with fitted 
absorbing column densities of (0.96±q j2) ^ W^^ cm~^ and 
(1.96±g-29) ^ ^^^^ cm~^ respectively. The metallicity for the two 
thermal components were fitted simultaneously giving a value of 
0.14 ± O.O1Z0. The fitted absorbing column density and photon 
index of the power-law component were (3.74±g ||) X 10^^ cm^^ 
and r = 1.43±[5;{]^ respectively. 

with the number of background sources that would be ex- 
pected within the D25 ellipse from 'The Chandra Deep Field 
South Survey' (Campana et al. 2001) for our completeness 
limit (~ 2 for our completeness limit and up to ~ 7 down 
to a flux level corresponding to that of our faintest detected 
source) . 

The three sources detected in all three energy bands 
(Sources 9, 21 and 26), along with sources 11 and 29, 
which were only detected in the medium band, are most 
likely X-ray binaries as none of these sources are coinci- 
dent with the positions of the reported supernova explosions 
in NGC 5253 in 1895 and 1972 (1895B: a = 13''39'"57.4% 
S = -31°38'06"; 1972E: a = 13''39™52.6" , S = -31°40'10", 
Caldwell & Phillips 1989) In fact, none of the detected 
sources tie up with the reported positions. Source 21 lies 
closest to the reported position of SN1895B (18 arcsec to 
the SW of its position) and source 11 lies closest to that for 
SN1972E (15 arcsec WSW of its position). 

3.3 Diffuse Emission 

A background and point source subtracted spectrum of 
the region containing the diffuse emission was extracted 
from the Chandra data only, as the spatial resolution of 
XMM-Newton was inadequate for resolving the point sources 
within the central region. This diffuse emission spectrum is 
shown in Fig.|5| The fitted model shown is an absorbed three 
component fit having two thermal components and a power- 
law component - ■wahsGAL{wabs{mekal) + wabs{mekal) + 
wabs{po)). The fit gives Nh = (0.96±!];;i) x 10^^ cm"^ 
(1.96±g:it) X 10^^ cm^2 and (3.74±g:tl) x 10^^ cm^^ for 
the fitted absorbing column densities of the cool, warm and 
power-law components respectively. The fitted temperatures 
are kT ~ 0.24 ± 0.01 keV for the soft component, and 



kT ^ 0.75 ± 0.05 keV for the medium component, while 
the fitted photon index for the power-law component is 
r = 1.43±o;o7. The metallicities for the two thermal com- 
ponents were fitted simultaneously, with 0.14 ± O.OIZq. 

The absorption corrected fiuxes (and percentage con- 
tributions) in the 0.37 — 6.0 keV energy band in each of 
these three components were: cool thermal: (1.84±o;i3) x 
10"" erg s"^ cm"^ (20%); warm thermal: (1.47±g:li) x 
10""ergs-i cm-2 (15%) and power-law: (6.48±E!:g) x 
10"^"^ erg s~^ cm~^ (65%). The errors shown here are the 
90% confidence levels for the normalisation value obtained 
from the fit. This result suggests that the diffuse emission 
consists of at least two components at different temperatures 
and many unresolved point sources. 

As discussed in Section 3.1, several of the point sources 
in the central region of NGC 5253 could be young super- 
bubbles and the subtraction of such sources could therefore 
remove a substantial fraction of the diffuse emission. Com- 
paring the fluxes above with those obtained from the analy- 
sis of the total emission gives the result that the point source 
subtraction has removed ~ 25% of the cool thermal emis- 
sion, ~ 60% of the warmer thermal emission and 35% of the 
emission attributed to the power-law component, but these 
flgures do not take account of the fact that the fltted values 
for the absorbing column densities have changed between 
the two flts. Re-fitting the diffuse emission with the values 
fitted for the total emission and only allowing the normal- 
isations for the individual components to be free gives re- 
ductions in the fiuxes that correspond to ~ 12% for the cool 
thermal emission, ~ 46% for the warmer thermal emission 
and ~ 14% for the emission from the power-law component. 
This does suggest that the cooler thermal component is af- 
fected less by the point source subtraction and so probably 
lies further away from the point sources and star-clusters 
than the warmer component. When sources 16, 17, 18 and 
20 are not subtracted from the total emission and this mod- 
ified emission is fitted as above with only the normalisations 
of each component as free parameters, the reduction in the 
fiuxes for each component are: cool thermal: ~ 10%; warm 
thermal: ~ 43%; power-law: ~ 13%. These figures suggest 
there is emission from all three components of the fit as- 
sociated with these four sources, with maybe slightly more 
of the warmer thermal component than the other two, and 
yet the temperatures measured for them during the spec- 
tral fitting was of the order of that of the cooler component. 
As their contribution is not a substantial fraction of the to- 
tal emission in all three components, there inclusion in the 
source subtraction does not effect the fluxes measured for 
the diffuse emission greatly. 

Other parameters of the two gas components have been 
calculated, from the spectral flts to the Chandra data for 
the diffuse emission and the XMM-Newton data for the to- 
tal emission. The results are shown in Table 7, including 
the effect of a fllling factor /, and within errors most of the 
calculated values agree between the two data sets. The im- 
plications of these flgures for the evolution of NGC 5253 will 
be discussed further in Section 4. The assumption of spher- 
ical symmetry and a filling factor of unity are both likely 
to be overestimates, resulting in the quoted figures being 
underestimates for rie and P and overestimates for M, Eth 
and tcooi (see Strickland & Stevens 2000 for a discussion of 
filling factors and their likely values in galactic winds). 
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Table 5. Hardness ratios for the 17 sources detected within the D25 ellipse. Soft band, s, 0.3 — 1.1 keV , medium band, m, 1.1 — 3.8 keV 
and hard band, h, 3.8 — 8.0 keV . Where no counts are shown, wavdetect failed to detect the object in that energy band. Column 1 gives 
the source numbers as shown on Fig. Columns 2-4 are the counts in the 3 different energy bands and columns 5 and 6 give the values 
of the hardness ratios calculated as detailed in the column headings. 



Source 


Counts in 


Counts in 


Counts in 


(m - s) 


(h - m) 




Soft Band 


Medium Band 


Hard Band 


(m + s) 


(h + m) 




(0.3 - 1.1 koV) 


(1.1 - 3.8 keV) 


(3.8 - 8.0 keV) 






9 


31.1 ± 5.7 


36.1 ±6.2 


- 


0.07 ±0.13 


- 


10 




6.1 ±2.6 


7.9 ±3.0 




0.13 ±0.29 


11 


7.3 ±2.8 


12.7 ±3.7 




0.27 ±0.24 




13 




7.6 ±2.8 








15 


16.0 ±4.2 


13.0 ±3.7 




-0.10 ±0.19 




16 


28.0 ±7.1 










17 


27.1 ± 7.0 










18/19 


43.4 ±8.9 


28.7 ±6.5 




-0.20 ±0.16 




20 


29.7 ±7.6 










21 


182.2 ± 14.2 


190.9 ± 14.1 


23.0 ±4.9 


0.02 ±0.05 


-0.79 ± 0.09 


22 


63.3 ± 9.2 


12.2 ±4.1 




-0.68 ± 0.16 




24 












25 


7.5 ±2.8 


7.9 ± 2.8 




0.03 ±0.26 




26 


23.0 ±4.9 


20.0 ±4.6 




-0.07 ±0.16 




27 




34.4 ±6.0 








29 


9.4 ±3.2 


9.7 ±3.4 




0.02 ±0.25 





Table 6. Hardness ratios for the 17 sources detected within the D25 ellipse. Soft band, s, 0.3 — 0.8 keV , medium band, m, 0.8 — 2.0 keV 
and hard band, h, 2.0 — 8.0 kcV . Details as per Table 5. These bands are used for comparison with data for NGC4449. 



Source 


Counts in 


Counts in 


Counts in 


(m - s) 


(h-m) 




Soft Band 


Medium Band 


Hard Band 


(m + s) 


(h + m) 




(0.3 - 0.8 koV) 


(0.8 - 2.0 keV) 


(2.0 - 8.0 keV ) 






9 


17.3 ±4.2 


36.6 ±6.2 


17.1 ±4.4 


0.36 ±0.15 


-0.36 ±0.15 


10 






14.3 ±4.4 






11 




10.5 ±3.3 








13 












15 


8.9 ± 3.2 


17.7 ±4.4 




0.33 ±0.22 




16 


17.4 ±5.1 


17.3 ±5.4 




0.00 ±0.22 




17 


24.0 ± 6.3 


22.5 ±6.8 




0.03 ±0.20 




18/19 


23.6 ± 6.2 


44.5 ±8.6 




0.31 ±0.16 




20 


16.1 ± 5.3 










21 


66.2 ± 9.0 


264.2 ± 16.7 


73.0 ±8.7 


0.60 ±0.07 


-0.57 ±0.06 


22 


39.6 ± 7.2 


39.9 ±7.3 




0.00 ±0.14 




24 












25 


5.6 ± 2.4 


7.8 ±2.8 




0.16 ±0.28 




26 


16.7 ±4.1 


22.8 ±4.9 


6.1 ± 2.6 


0.15 ±0.16 


-0.58 ±0.22 


27 




17.0 ±4.2 


21.7 ±4.8 




0.12 ±0.17 


29 




12.9 ±3.6 









Overall, the metallicity of the hot X-ray emitting gas 
is similar to that obtained from optical spectroscopy of the 
central HII regions (Webster & Smith 1983) of the galaxy. 
Compared to other dwarf galaxies though, the values of 
O.14±O.O1Z0 from the diffuse X-ray data and 0.19±Ei:[]i Zq 
from the optical data are low. Typical values for some 
other dwarfs are: Mrk 33: ~ 0.3^0 (Legrand et al. 1997); 
NGC4449: ~ O.32±O.O8Z0 for the X-ray emitting gas (Sum- 
mers et al. 2003); NGC 1569: > 0.25^0 for the X-ray emit- 
ting gas and a mean value of 0.2Zq for the ISM material 
(Martin et al. 2002). In a further attempt to investigate 
the metallicities of the thermal components, and in partic- 
ular the ratio of the a-elements to Fe, the spectrum was 



re-fitted using variable abundance mekal models {vmekal) 
for the two thermal components. The abundances of the 
two thermal components were tied and the individual abun- 
dances of Mg, Ne, Si and Ca relative to Solar values were 
tied to that of O to form the group of a-elements. The abun- 
dances of all the other elements were set to the global value 
determined above. The a/Fe ratio obtained in this way was 
found to have no significant difference to the Solar value 
at (1.05 ± O.3O)(a/7^e)0. The confidence contour plot for 
these fits is shown in Fig. 1101 showing that the a and Fe 
abundances are highly correlated, and that the fitted ratio 
is ~ 1. We note that this is somewhat in contrast to the val- 
ues of Martin et al. (2002) for the dwarf starburst NGC 1569, 
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Table 7. The gas parameters for the two thermal components of the diffuse emission. The assumptions made are that the volume of the 
emitting region is V = 1.54 X 10®* cm^ (assuming spherical symmetry); D = 3.15 Mpc; a filling factor, /; rie ~ (EI /V f)^^'^ where EI is 
the emission integral (normx47rZ)^)/10~^* and norm is the normalisation obtained from the spectral fitting; P ~ 2nekT, M ^ riempVf, 
^Th ~ 3nf,kTV and tcoo! ~ (3fcT)/(Ane) where A = Lx/EI. Columns 2 and 3 give the values calculated from the Chandra data while 
columns 4 and 5 are the corresponding XMM-Newton values. 



Chandra XMM-Newton 





Soft component 


Medium component 


Soft component 


Medium component 


kT ( keV ) 


0.24±0.0i 


0.75±0:«| 




fi7+0.04 
"■"'^0.05 


T (lO^K) 




8-69±E5:i 






Lx (lO^S erg s"! ) 


2 1 q+0.16 
^■1^^0.15 


1 7>;+0.l5 


1 fic;+0.08 


1.82±0:14 


ne [X/V2] (cm-3) 


U.UOOztg Q()2 


n 041+0.002 


ofio+o-foi 


044+0'002 


Eth [x/1/2] (1054 erg) 




o 9K_l_0.18 




2.2i±g:li 


M [x/i/2] (io6 Mq ) 


0.84±0.03 


KQ+0.03 






P [xj-i/^] (10-11 dyn cm-2) 


4 Qv+o.ao 


n 7Q+0.79 


4.82±0:i0 




tcool [x/1/2] (lOSyr) 


1-°'='=0.19 


*-l'^='=0.56 


2 ic;+0.16 


3.82±g:| 


A (10-24 erg s-l cm^) 


3.36±q;33 


6.78±{]-?| 


2.95±g;22 


5.95±Q;7g 


EI (lO'^lcm-3) 


a cn_i_0.45 


2.58±o;i 


5.59±0:i 


3.06±0:i 



where values of the a to Fe abundance ratio of 2.1 — 3.9 times 
the Solar value are quoted (and in the case of NGC 1569 
the photon statistics are rather better than for NGC 5253). 
This value for NGC 5253 is more in line with the value of 
(0.91±[l:?;^)(a/Fe)Q found for the dwarf starburst NGC 4449 
(Summers et al. 2003) and may be indicative of differences 
in the contributions made to the hot gas from SNe in these 
different starbursts. The high values quoted for NGC 1569 
suggest the presence of more Type II SNe, which would be 
expected if its current starburst phase is ~ 10 — 20 Myr old 
as suggested in the literature on this starburst (e.g., Israel 
& De Bruyn 1988; Gonzalez-Delgado et al. 1997; Hunter et 
al. 2000). Such activity would give rise to more a-elements, 
as would a pattern of shorter more intense bursts of star 
formation. The fact that NGC 5253 has shown a fairly uni- 
form rate of star-formation throughout the galaxy over the 
last 10* - 10^ yr (CaldweU & Phillips 1989), would suggest 
that, with the exception of the current burst in the very cen- 
tral region, which could be very young and hence have little 
Type II SNe activity at present, its star- formation and en- 
ergy injection from star-forming regions has been generally 
of a more continuous nature. 
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Figure 10. The confidence contours for the cx and Fe element 
abundances (see text for details). The units of both axes are the 
solar values for the a elements and Fe respectively. The contours 
shown are the 68%, 90% and 99% contours. 



From Table 4, it can be seen that there is little tem- 
perature variation evident within the sources, that appear 
to be extended (16, 17, 18, 19 and 20), in the central region 
of NGC 5253, with the exception of source 16 in the SW 
being cooler than the rest. The absorbing column densities 
fitted to them are also fairly consistent and there is overlap 
seen between adjacent sources within errors. The dust lane 
in NGC 5253 bisects the nuclear region such that sources 18 
and 19 lie to the north of it and sources 16, 17 and 20 are to 
the south. From the hardness ratios and source counts in var- 
ious energy bands shown in Tables 5 and 6, the two northern 
sources are harder than those in the south, suggesting that 
they should be either hotter or lie behind larger absorbing 
columns. Source 17's temperature and 16's absorbing col- 
umn density would argue against this. The latter could be 
explained by the western side of the galaxy being tilted away 
from us. The low count-rate statistics however, hamper the 
drawing of any firm conclusions, though the results do sug- 
gest that these objects are similar in nature to each other. 



Also, from the 3 colour Chandra image, further differ- 
ences can be seen in the nature of the point sources in these 
three energy bands. Sources 17, 18, and 20 all appear very 
red on the three colour image, reinforcing their similar na- 
ture. Sources 16 and 19 show slightly different behaviour, by 
having more emission in the medium band than the other 
three sources, which is reflected in their yellow colour on 
Fig. 0] The diffuse emission is mainly confined to the cen- 
tral region of the galaxy, and the peak of its emission is 
centred on the region occupied by sources 17, 18 and 19, 
which also corresponds to the region of peak intensity in the 
Ha emission. The images of the XMM-Newton data shown 
in Figs. |5| and |3 show that the diffuse X-ray emission ex- 
tends away from the central region in several directions, al- 
though as stated earlier most of the extension to the SW 
appears to be associated with point sources. This is con- 
firmed by the log-log plot of the radial surface brightness 
profile shown in Fig. 1111 produced by extracting the counts 
from 20 concentric annuli each with a width of 2.75" cen- 
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tred on a position coincident with the centre of the D25 
ellipse at Q = 13''39™56.2'' and S = -31°38'29.9". This 
figure shows that the difluse emission can only be traced 
out to a distance of ~ 35" (0.53 kpc for our assumed dis- 
tance) before it falls to a level comparable to that of ran- 
dom background fiuctuations. The preferred directions for 
the observed extensions in the diffuse emission seem to be 
along the minor axis of the galaxy, both to the NW and 
SE. If these extension in the diffuse X-ray emission are due 
to superbubbles and superwinds then their activity should 
have swept-up shells and produced filaments which should 
be evident in the Ha image of the galaxy. To look for such 
evidence, of correlations between the extended X-ray emis- 
sion and structure within the Ha emission. X-ray contours 
from the point source and background subtracted diffuse 
emission have been overlaid on the Ha image of the central 
region of NGC 5253. In addition, the detected sources in the 
central region of the galaxy (sources 15-22) are also over- 
laid on the image to show the relative positions, and extent 
of these sources (as determined by the wavdetect program), 
compared with the regions of increased star-formation on 
the Ha image. The results are shown in Fig. 1121 and the 
correlations between the various components are discussed 
further in Section 4.1.1. As a further comparison between 
emission at different wavebands, Fig. ll3l shows the positions 
of the six largest star clusters found in the central region 
of NGC 5253 overlaid on the Chandra soft (0.3 - 1.1 keV) 
and medium (1.1 — 3.8 keV ) smoothed images of the galaxy 
along with the positions of the X-ray sources 15-22 as shown 
in Fig. 1121 These images show only two strong correlations 
between the most northern and southern star clusters shown 
here and X-ray sources. Most of the soft X-ray emission oc- 
curs between and around these large star clusters in regions 
where many smaller clusters are to be found. The second 
most northern of these large star clusters shown on these 
images has four smaller clusters in close proximity to it, yet 
it appears to have little X-ray emission close to it. However, 
these clusters are located at the western end of the dust 
lane, which is most likely responsible for the absorption of 
the X-ray emission seen both in these images and on Fig. |1| 



4 DISCUSSION 

It is apparent from Fig. I12l that both the diffuse X-ray emis- 
sion and the Ha emission are mainly confined to a region 
of radius ~ 0.5 kpc in the centre of NGC 5253. The most 
intense X-ray emission is coincident with the most intense 
emission seen in the Ha image in the Northern half of the 
nucleus, suggesting that this emission is associated with the 
increased star-formation occurring in these regions. In gen- 
eral, higher column densities were fitted for the sources to 
the W and S of the galaxy, as seen in Table 4, and would sug- 
gest that these areas lie further away from us behind more 
absorbing material than the N and E of the galaxy. This 
would mean the minor axis lies away from our line-of-sight 
(tilted both to the N and W), and that the galaxy is also 
rotated about this axis in such a way that the Southern side 
of the galaxy disk is tilted away from us, consistent with 
NGC 5253 being a dwarf elliptical galaxy that is inclined to 
our line-of-sight. The bulging of the diffuse emission seen to 
the NW and SE in the X-ray contours and XMM-Newton 
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Figure 11. Background subtracted radial surface brightness plot, 
showing the extent of the difEuse X-ray emission in the 0.3 — 
8.0 keV band. The counts were extracted in 20 concentric annuli, 
each with a width of 2.75" centred on a = 13''39'"56.2'' and 
S = -31°38'29.9". The diffuse emission extends out to ~ 35" 
(~ 0.53 kpc) by which point, the emission has dropped to the 
background level. 



images, could be indicative of a bipolar outflow along the 
minor axis of the galaxy, as would be expected from stan- 
dard superbubble models (Weaver et al. 1977). 



4.1 Comparison with Observations at Other 
Wavelengths 

4.1.1 Ha Emission 

NGC 5253 contains many regions of increased star-formation 
within the central region of the galaxy, which are clearly vis- 
ible in Ha observations, with the most intense regions being 
seen in the North of the nucleus, coincident with the end of 
the dust lane and weak CO emission (Turner et al. 1997). 
This is the region that contains two of the brightest star- 
clusters analysed by Tremonti et al. (2001) and correspond- 
ing to our detected X-ray sources 18 and 19. The extension 
seen in the diffuse emission to the WNW of this region could 
be the result of an outflow induced by the stellar winds of 
the massive stars present in these star clusters. This region 
also lies close to the position of the 'radio super-nebula' dis- 
cussed by Gorjian et al. (2001), which they have interpreted 
as an highly obscured, very young globular cluster, rather 
than an obscured AGN (see also Mohan, Anantharamaiah 
& Goss 2001; Turner et al. 2003). 

As well as many of the star clusters coinciding with X- 
ray sources, some of the filamentary and shell- like structures 
seen in Ha images seem to have counterparts in the diffuse 
X-ray emission. The small Ha shell seen to the SSE in Fig. 1121 
appears to be filled with diffuse X-ray emission and the ra- 
dial filament extending off the Ha image to the ESE is also 
associated with what appears to be a filamentary extension 
of the diffuse X-ray emission. The field-of-view of the Ha im- 
age shown in Fig. I12l is too small to show the ~ 1 kpc shells 
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Figure 12. X-ray contours from the Chandra diffuse emission (0.3 — 8.0 keV energy band) overlaid on an Ha image of NGC 5253. The X- 
ray contours shown are at flux densities of 0.26 X lO"'^^, 0.52 X 10~^^, 1.04 X lO"'^^, 2.08 X 10~^^ and 4.16 X lO"'^^ erg s~^ cm~-^ arcmin"-^. 
Also shown are the positions and extent of the sources detected in the central region of the Chandra data and these are labelled with 
their ID numbers as shown in Fig. Q North is to the top and East to the left. 



detected to the West of the galaixy by Marlowe et al. (1995), 
but these can be seen in Fig. 2j of Martin (1998). The ori- 
entation of the more northern of these two shells is roughly 
aligned with the galaxies minor axis. The slow expansion 
speed of these shells [~ 35 km s"'^ (Marlowe et al. 1995)] 
and the fact that the diffuse X-ray emission does not extend 
out to them, as shown in the radial surface brightness plot 
of Fig. 111! may mean that these shells are not associated 
with the current burst of star formation but are remnants 
of earlier starburst activity. The presence of more than one 
bubble having been blown close to the minor-axis on the 
West side suggests that what is being seen in this galaxy 
is much more complex than the simple bipolar outflow sce- 
nario assumed to explain the effect of starburst activity. This 
phenomenon could be the result from a young starburst, in 
so much as the outflows from individual star-clusters may 
not have overlapped or merged to produce what would be 
observed as a single outflow. The presence of the dust lane 
bisecting the nuclear region may also be acting in such a way 
as to channel any superbubbles so that they expand prefer- 
entially into areas where the density of the ambient ISM 



is less. The possible dynamics of such superbubbles will be 
discussed later. 



4-1.2 Infrared and Radio Emission 

Strong mid-infrared emission (Gorjian et al. 2001) is de- 
tected from a source in the northern half of the central re- 
gion of the galaxy just at the Northern edge of our detected 
sources 18 and 19. This is an extremely compact source oc- 
cupying a region with diameter ~ 1 — 2 pc and needing to 
contain ~ 10^ — 10^ stars to supply its ionization require- 
ments. This object is most likely a very young (no more than 
a few X 10^ years old) globular cluster. 

Just to the south of this infrared emission, at a position 
corresponding to our source 19, is the peak of the 2cm radio 
emission (Beck et al. 1996). While mainly concentrated in 
the Northern half of the central region of the galaxy, some 
2cm emission is also detected at the positions of our sources 
17 and 20 but not coincident with source 16. This radio emis- 
sion has a flat spectrum, with a spectral index of — O.liO.Ol, 
which is attributed to thermal bremsstrahlung emission from 
the HII in the region. Relatively low levels of non-thermal 
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synchrotron emission are seen in NGC 5253 (Turner et al. 

1998) , which is unusual for a starburst galaxy and again 
suggests that the starburst occurring in the central region 
of NGC 5253 could be young. 

As mentioned earlier, VLA 21 cm observations of 
NGC 5253 have suggested that the bulk of the HI rotates 
about the galaxy's major axis (though see Meier et al. 2002 
for evidence of infall into NGC 5253) . Fig. 1 of Kobulnicky 
& Skillman (1995), shows that the neutral gas extends to 
about twice the optical extent of the galaxy along its mi- 
nor axis, occupying a region of ~ 4' in diameter. This is 
well beyond the extent of the diffuse X-ray emission and Hq 
emission seen in the central region. The peculiar rotation 
could result from large-scale outflow along the minor axis, 
but such behaviour should give rise to a shell-like centrally 
evacuated structure in HI emission, as seen in the ~ 1 kpc 
Ha shells. The column density measured for the HI increases 
towards the centre and peaks at a position close to our de- 
tected source 20. This is offset by about 15" from the peaks 
of emission seen at shorter wavelengths and is in the south- 
ern rather than northern half of the galaxy, however this 
does coincide with the point where the dust lane enters the 
nuclear region from the ESE. 

4.2 Potential for Blow-out and its Implications on 
the Fate of NGC 5253 

Although at present there does not appear to be a clearly 
defined single superbubble structure within NGC 5253, the 
presence of hot, over-pressured gas with a long cooling time, 
should lead to an outflow from the central region. With this 
in mind, in this section we apply the standard superbubble 
model to the hot gas of NGC 5253 to investigate the poten- 
tial of the galaxy to deposit energy, metal-enriched material 
and ISM material into the intergalactic medium. 

The age of the current burst of star-formation in the 
centre of NGC 5253 is not well determined, and probably 
does not have a single age. The comparative lack of a non- 
thermal signature in the radio-spectrum (Beck et al. 1996; 
Turner et al. 1998) may suggest a lack of supernova activ- 
ity and the detected presence of Wolf-Rayet stars (Schaerer 
et al. 1997) in regions corresponding to our detected source 
19 suggest that star-formation has been occurring for only 
a few Myr in this region. Typically ages for star clusters 
in NGC 5253 reported in the literature are: 7.9 Myr (Mar- 
lowe et al. 1999); 2.8 Myr and 4.4 Myr for two separate 
regions within the nucleus (Schaerer et al. 1997); 2 Myr 
for the brightest star cluster (Tremonti et al. 2001); 5 Myr 
(Martin & Kennicutt 1995); 3 Myr (Mas-Hesse & Kunth 

1999) , though as noted earlier there are some older clusters 
(Calzetti et al. 1997), and indeed these older star-clusters 
can still be contributing mechanical luminosity. 

Based on these figures we adopt an age of 5 Myr for the 
current starburst, in order to perform some analysis of the 
dynamics of the probable superbubbles in NGC 5253. We 
recognise the limitation of this approach in such a complex 
system as NGC 5253, where the star-formation is not co- 
eval. The northern star-clusters in NGC 5253 are younger 
and here the analysis may be useful. 

The bulk of the diffuse X-ray emission is confined within 
a radius of ~ 0.6' (0.53 kpc for our assumed distance), that 
corresponds to the same region containing the majority of 



the Hof emission as shown in Fig. 1121 For simplicity this will 
be treated as a single superbubble resulting from the com- 
bined action of the multiple star-clusters within the star- 
burst. The dynamic time required for the two ~ 1 kpc Ha 
shells detected by Marlowe et al. (1995) to have evolved is 
~ 10 Myr, assuming they have had a constant expansion 
velocity of 35 km s~^ (Marlowe et al. 1995) for this time. 
This value is greater than the estimated duration of the cur- 
rent burst of star-formation, and the diffuse emission does 
not appear to extend this far out. However, this figure is 
also likely to be an over estimate since their speed will have 
been greater in the past, which would bring their expansion 
times closer to the higher age estimates of the star clusters. 
These observations raise the question of what is responsible 
for their existence which will be discussed further later. 

The combined action of the stellar winds of the massive 
O and B stars in the star-forming regions in the centre of 
NGC 5253 will blow a superbubble that can be modelled us- 
ing the starburst driven outfiow models of Castor, McCray 
& Weaver (1975) and Weaver et al. (1977). In such mod- 
els, the stellar winds from the OB associations sweep up 
and shock a thin, dense Ha emitting shell of ISM material, 
with a contact discontinuity just behind the shell, separat- 
ing the shocked wind material from the shocked ISM ma- 
terial. The innermost region of the bubble will contain the 
freely expanding supersonic wind bounded by an inward- 
facing shock. In the 'snow-plough' phase, the shocked wind 
material occupies most of the bubbles interior and most of 
the bubble's mass is contained in the Ha shell. The expan- 
sion of the bubble is driven by the over-pressured interior, 
which, as can be seen from the figures of Table 7 is of an or- 
der of magnitude greater than typical values for the ambient 
ISM within the Milky Way which is ~ 10"^^ dyn cm"^ The 
radius Rb and expansion velocity fs, of the bubble during 
this phase, are given by 

/ T \ 1/5 

Rb = 168 tJ pc (1) 

V no J 

/ r \ 1/5 

nn / ^mec/i. 40 \ ,—2/5 i —1 /n\ 

«s = 99 tg ' km s , (2) 

V no / 

where LmechAO is the mechanical luminosity injected by the 
starburst in units of 10*" erg s~^ , ta is the age of the star- 
burst in units of 10^ yr and no is the ambient atomic number 
density of the ISM. 

The mechanical luminosity (Lmech) injected by the star- 
burst can be estimated from the stellar populations in the 
starburst. Marlowe et al. (1999) have performed a detailed 
analysis of the starburst in NGC 5253, using starburst evo- 
lutionary synthesis models (c.f. Leitherer & Heckman 1995; 
Leitherer et al. 1999) and they estimate that the energy 
injection rate in NGC 5253 is Lmech ~ 3.6 x 10*° erg s"'^ 
for both instantaneous burst and continuous star-formation 
models, and we shall adopt this figure hereafter. The mass 
injection rate from the Starburst99 models for this level of 
mechanical luminosity will be ~ 0.02 — 0.04 Mq yr~^ de- 
pending on whether instantaneous or continuous starburst 
models are assumed (for an age of 5 Myr). 

An estimate of the ambient density in the central re- 
gions, no, can be made based on the HI observations of 
Kobulnicky & Skillman (1995). The peak HI column den- 
sity in the central regions is 2.6 x 10^^ cm~^ which, assum- 
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Figure 13. Positions of the six largest optical star clusters are shown as crosses and these are overlaid on the Chandra soft-band 
(0.3 — 1.1 keV ) smoothed X-ray images of NGC 5253. Also overlaid on these images are the X-ray sources detected in the central region 
of the Chandra data, which are labelled in Fig. 1121 N is to the top and E to the left and the image is a linear greyscale. The flux densities 
range from to 7.0 X 10~^^ erg cm~^ arcmin"^. 



ing that the path length is ~ Ikpc (i.e comparable to the 
beam width) leads to no = 1 cm~^. We note that Martin 
(1997) inferred densities of ~ 250 cm~^ in the warm pho- 
toionised gas in the core of NGC 5253, though with a filling 
factor of a percent or so, and that much higher densities are 
also implied for the smaller volumes associated with the HII 
regions in NGC 5253 (Turner et al. 1998). 

Using these figures in equations 1 and 2 give values for 
Rb and vb of 570 pc and 67 km s~^ . The size estimate is 
comparable to the size of the region occupied by the bulk of 
the X-ray and Ha emitting gases, but no allowance has been 
made for the effect of the interaction of outflows from the 
multiple star-clusters present in the central region, while the 
theory is applicable to the growth of a single bubble. This 
size estimate is smaller than the size of the reported large- 
scale Ha shells, their expansion velocities were measured 
to be ~ 35 km s~^ , which is of roughly the same order as 
predicted here. It is possible that we may be observing the 
escape of some stellar wind material through ruptures in 
the original bubble that are then proceeding to sweep-up 
and shock more of the ISM lying outside the central re- 
gion. The emission from such gas has to be at a very low 
intensity though as it is not detected by either Chandra or 
XMM-Newton. The average count rate in the annuli beyond 
the 0.53 kpc to which the X-ray emission can be traced is 
0.05 ± 0.28 X 10""^ cts s~^, which corresponds to a flux of 
~ 0.3 ± 1.9 X 10~^^ erg s~^ cm~^. Alternatively, a clumpy 
environment could allow preferential escape along channels 
between molecular clouds, although there does not appear 
to be evidence for the existence of such objects, and the 
low expansion velocities observed would argue against such 
uninhibited outflow. There are however two 'dark clouds' 
reported in the West of the galaxy by Hunter (1982). 

The superbubble model allows prediction of the X-ray 
luminosity of an expanding bubble to be made. From Ta- 



ble 7, it can be seen that the radiative cooling time for both 
components of the X-ray emitting gas is long compared to 
our assumed expansion time of 5 Myr, and so radiative losses 
can be assumed to be small. Assuming spherical symmetry 
for the expansion, uniform density for the ambient ISM and 
a constant rate of kinetic energy injection, the X-ray lumi- 
nosity is given by 

Lx^ J n{rfAx{T,Z)dV (3) 

where Ax {T, Z) is the volume emissivity of the gas be- 
ing considered. Based on the density n(r) and temperature 
T{r) profiles (Chu & Mac Low 1990), this equation can be 
evaluated to yield the X-ray luminosity of the superbub- 
ble. Assuming an emissivity of the X-ray emitting gas of 
A = 6.8 X 10"^'' erg s"^ cm^ (cf Table 7), we find that for 
NGC 5253 

r no, ir,36r33/35 17/35,19/35 -1 

Lx = 9.3 X 10 L„i,;,,4ono ' ' erg s , (4) 

As discussed earlier, assuming an energy injection 
rate from the starburst in NGC 5253 of Lmech = 3.6 x 
lO^'' erg s~^ , an ambient density of no = 1 cm"'' and an 
age of 5 Myr, leads to a predicted thermal X-ray lumi- 
nosity for NGC 5253 of Lx = 0.7 x 10^* erg s'^ , which is 
lower than that observed for the two thermal components 
{Lx ~ 4 X 10'^* erg s~^ , but given the uncertainties in the 
some of the assumed parameters, and indeed the validity of 
treating the situation in NGC 5253 as a single superbubble, 
this is perhaps reasonable agreement. 

We note that an estimate of the mechanical energy in- 
jection rate can also be obtained from the thermal energy 
contained within the diffuse X-ray emission. Averaging the 
values obtained from the Chandra and XMM-Newton data, 
the total thermal energy is 3.4 x 10^* erg. Assuming this is 
the result of complete thermalization of the stellar winds 
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for the duration of the current starburst, then Lmech is 
^ 2.2 X 10'"^ergs~^, which is again in broad agreement 
with that from theoretical considerations of the starburst 
stellar population. 

The total mass of material contained in the X-ray emit- 
ting gas obtained from the average figures of the two data 
sets is ~ 1.4 X 10® M© (ignoring the filUng factor). If this 
mass has been injected during the lifetime of the starburst 
then the mass injection rate is on average ~ 0.3 M© yr~^ . 

Another estimate of the mass injection rate can be ob- 
tained from the fact that for an adiabatic superbubble, the 
X-ray temperature Tx ~ 5/llTo, where 

2/L_^X/im^ ^ 
3 \ M J k 

(Heckman et al. 1995) and is the temperature attained if all 
the mechanical energy injected by the starburst remains en- 
tirely within the shocked wind region. The weighted average 
temperature of the diffuse X-ray emitting gas (Table 7) is 
~ 5 X 10" K, which gives a value for To ~ 1.1 x 10^ K 
and from the assumed Lmech, a mass injection rate of 
M ^ 0.25 M© yr~^ . Finally, an estimate of the deposited 
mass from both the stellar winds and conductive evapora- 
tion of the swept-up ISM material can be made from the 
density of the difi'use X-ray emitting gas. This figure de- 
pends on the filling factor of the bubble interior but assum- 
ing / = 1 results in an upper limit. Further, it is necessary 
to assume that any volume occupied by the swept-up cool 
shell of ISM and the freely expanding stellar winds is small in 
comparison to the shocked wind region. With the combined 
density of the soft and medium components, n ~ 0.1 cm~^, 
V = 1.54 X 10®* cm^ and all other symbols having their 
usual meanings, 

M = pV = nuniH^nR^ ~ 8 x 10® M© , (6) 

o 

which implies a mass injection rate of 0.16 M© yr~^ over 
5 Myr. 

The mass injection rates into the hot phase derived here 
from the X-ray observations are thus either slightly above 
or broadly comparable to the current star-formation rate 
(0.2 M© yr^^ ) derived from infrared and Hq luminosities 
(Table 1), but rather higher than the current stellar mass- 
injection rates. 

The mass deposition rate into the hot-phase seen here 
suggests that if the outflow from the starburst region de- 
velops into a supcrwind and escapes the gravitational po- 
tential well of the galaxy then there is substantial scope for 
both metal-enrichment of, and deposition of large amounts 
of mass and energy into the IGM. The likelihood of blow-out 
occurring can be assessed by applying the criterion for blow- 
out, defined by Mac Low & McCray (1988) to NGC 5253. 
This defines a parameter A, the dimensionless rate of ki- 
netic energy injection, and predicts blow-out if A > 100, 
with A given by 

A = 10^L^ech,40Hfc"pe^'4"''"'^J'" (7) 

where Lmech, 40 is the mechanical energy luminosity in units 
of 10*" erg s^^ , Hkpc is the galaxy scale-height in kpc, and 
P4 is the initial pressure of the ISM in units of P/k = 
10* K cm"*. For A > 100, with Lmech ao = 3-6, no = 5 cm"* 
and assuming P4 ~ 1 (typical of the value in the Milky 



Way), then Hkpc < 9. The HI halo of NGC 5253 extends to 
around twice the optical extent ~ 2' 1.83 kpc) along the 
minor axis of the galaxy (Kobulnicky & Skillman 1995). It 
seems extremely likely then that the hot, metal-enriched gas 
generated by the current starburst along with the swept-up 
and conductively evaporated ISM material can escape the 
gravitational potential well of NGC 5253. 



5 SUMMARY AND CONCLUSIONS 

In summary, we have presented an X-ray analysis of the im- 
portant dwarf elliptical, starburst galaxy NGC 5253 using 
data from both the Chandra and XMM-Newton satellites. 
We detected X-ray emission from 31 discrete point sources 
in the Chandra ACIS-S3 chip data with 17 of these lying 
within the optical extent of the galaxy, as measured by the 
D25 ellipse. Some of these can be clearly identified as X-ray 
binaries while others seem to be associated with the emis- 
sion from super star-clusters rather than individual point 
sources, as they are extended compared to the Chandra PSF. 
The 5 sources in this category (sources 16 — 20) all lie in the 
central region of the galaxy and are associated (particularly 
in the case of sources 17, 18 and 19) with the most intense 
regions of emission seen on the Ha image. The unabsorbed 
luminosities of these sources in the 0.3 — 8.0 keV energy 
band, range from (0.25 — 1.50) x 10'^*^ erg s"'^ . 

The majority of the X-ray emission is confined to the 
central region of the galaxy where the intense, current burst 
of star-formation is occurring. The extent of the emission 
also correlates well with that of the Hq emission and both 
are confined to a region of diameter ~ 1 kpc in the centre 
of the galaxy. Both of these emissions are embedded within 
a spherical distribution of HI that extends out to ~ 2 kpc 
from the centre. 

The diffuse X-ray emission in NGC 5253 requires two 
thermal and one power-law component to best fit its spec- 
trum, suggesting that the hot gas is a multi-phase environ- 
ment with unresolved point sources, probably lower luminos- 
ity X-ray binaries, embedded in it. The fitted gas tempera- 
tures are 0.24 ± 0.01 kcV and 0.75 ± 0.05 keV for the soft 
and medium components respectively and their respective 
absorption corrected luminosities in the 0.37 — 6.0 keV band 
are (2.19±8:ii) x 10*® erg s"^ and (1.75±g:}t) x 10*® erg s"^ . 
This emission is seen to be extended in several directions in 
the Chandra data but this is not confirmed by the XMM- 
Newton data. However, the extent is not as great as the 
~ 1 kpc shells reported to the NW and WSW by Marlowe 
et al. (1995), in the Ha emission. These shells are expand- 
ing with velocities of ~ 35 km s"^ and hence have required 
~10 Myr to have evolved. NGC 5253 is proving to be a rather 
enigmatic object though, as our analysis of the dynamics of 
the hot gas are not consistent with the expansion velocity 
and age of these Ha shells. 

The hot X-ray emitting gas has a total thermal energy 
content of 3.4 x 10** erg and a total mass of 1.4 x 10® M© . 
From values in the literature, we assume an age for the star- 
burst responsible for this hot gas of 5 Myr and a value for 
the density of the ambient ISM of no — 5 cm"*. The ra- 
dius of the region in which most of the diffuse X-ray and 
Ha emission is confined is 0.53 kpc. Application of stan- 
dard superbubble models to NGC 5253's starburst, assum- 
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ing the age and density above, result in a predicted radius 
and expansion velocity for a single superbubble of ~ 400 pc 
and ~ 50 km s"'^ , figures which are comparable to those 
observed. 

The main problem hero would seem to lie in the fact 
that no clear bipolar outflow is observed in NGC 5253 and 
no allowance has been made for the presence of multiple 
superbubbles in the central region. If the central staxburst 
is young, the bubbles blown around individual star clusters 
will overlap and impede each others expansion confining the 
emission to a smaller region than expected and reducing the 
observed expansion rate. The Ha shells, reported by Mar- 
lowe ct al. (1995), to the West could be the result of the 
superbubbles closest to the edge of the central region being 
able to expand more freely. The lower velocity seen could re- 
sult from less stars contributing to the energy input in these 
regions and/or the outflow having experienced some recent 
deceleration on encountering higher density material such as 
the dark clouds reported in the West of the galaxy (Hunter 
1982). However, our radial surface brightness profile does 
not show any detectable X-ray emission beyond ~ 0.5 kpc. 

NGC 5253 is a very complex object, which shows that 
standard superbubble models are difficult to apply to dwarf 
starbursts galaxies. However, the energy injection rate into 
the galaxy would seem sufficient to allow the expanding hot 
gas to escape the gravitational potential well of NGC 5253 
and its relatively small (~ 1.8 kpc along the minor axis) 
HI halo would also offer little resistance to this. Tims, it 
seems likely then that NGC 5253 will lose metal-enriched 
material, mass and energy as a result of its current bout of 
star-formation. 
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